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We present a comprehensive first-principles study of the electronic charge redistribution in 
atomically sharp LaAlO3/SrTiO3(00f ) heterointerfaces of both n- and p-types allowing for non- 
stoichiometric composition. Using two different computational methods within the framework of the 
density functional theory (linear combination of atomic orbitals and plane waves) we demonstrate 
that conducting properties of LaAlO3/SrTiO3(001) heterointerfaces strongly depend on termina- 
tion of LaAlOsfdOl) surface. We argue that both the "polar catastrophe" and the polar distortion 
scenarios may be realized depending on the interface stoichiometry. Our calculations predict that 
heterointerfaces with a non-stoichiometric film — either LaO-terminated n-type or A102-terminated 
p-type — may exhibit the conductivity of n- or p-type, respectively, independently of LaAlC>3(001) 
film thickness. 
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I. INTRODUCTION 

The discovery of conducting interfaces between two ini- 
tially insulating materials — Ti02-terminated (001) sur- 
face of SrTi0 3 (STO) substrate and LaA10 3 (LAO) thin 
film deposited on top of it 1 — has attracted strong sci- 
entific interest during the last few years. 2-6 The high 
application potential of LAO/STO heterointerfaces has 
been demonstrated, e.g., by fabrication of highly voltage- 
tunable oxide diodes 7 that utilize the advantage of the 
electric-field controlled interfacial metal-insulator tran- 
sition of LAO/STO. 8 ' 9 

Conductivity of an atomically flat interface in the limit 
of large film thickness can be understood 1 from electro- 
static considerations within the so-called "polar catas- 
trophe" picture. From the perspective of formal charges, 
the atomic planes in [001] direction (which we refer to as 
monolayers) are neutral for STO (SrO° and TiO°), but 
charged for LAO (LaO + and A10^~). Transition from 
STO to LAO can be p-type (from SrO° to AlOj) or 
n-type (from TiO° to LaO + ). The corresponding jump 
of the surface charge at the interface would create an 
electric field inside LAO increasing linearly with the dis- 
tance from the interface — a "polar catastrophe". This 
bulk polarization of the LAO film can be compensated 
(thus averting the "catastrophe") if 0.5e per unit cell 
area is transferred from the LAO film surface onto the 
interface, resulting in a maximal sheet carrier density of 
n = 0.5/a 2 = 3.3 x 10 14 cm~ 2 (here a = 3.90 A is the lat- 
tice constant of STO assuming epitaxial matching of the 
LAO film). This estimate is immune to dielectric relax- 
ation and bond covalency/charge smearing effects 6 and 
thus provides a useful reference in the thick film limit. 

For sufficiently thin films, however, the polar catas- 
trophe may be tolerated 10 and a metal-insulator tran- 
sition occurs 11 as a function of the number of epitax- 
ial monolayers of LAO deposited. From electrostatic 
perspective, in a sufficiently thin film the internal field 



does not develop a potential difference large enough to 
overcome the dielectric gap. The accumulation of op- 
positely charged monolayers leads to progressive band 
bending until the critical thickness (5 u.c. or 10 mono- 
layers for n-type structures, according to experimental 11 
and theoretical 12-14 evidence) is reached beyond which 
the chemical forces are overcome, and the charge redis- 
tribution occurs. This mechanism is known as "polar 
distortion". 4 ' 10 As far as it is known, p-type interfaces 
do not exhibit this mechanism, as covalent forces over- 
whelm electrostatic ones. 

Similar electostatic arguments may be applied to non- 
stoichiometric LAO/STO structures, i.e. the ones with 
an odd number of LAO monolayers. The LAO films in 
these structures possess one extra electron (or a hole for 
p-type interfaces) per unit cell area compared to the par- 
ent bulk material, thus they should be conducting (with 
n = I /a 2 = 6.6 x 10 14 cm -2 ) irrespective of the thickness 
or the presence of STO substrate. However, as seen from 
the stoichiometric example, for thin films the competition 
between the semi-covalent bonds and long-range electro- 
statics is very sensitive to the the number of monolayers 
deposited. The nature of the conducting layer in non- 
stoichiometric LAO/STO interfaces is the main subject 
of our ab initio investigation. 

Recent experimental reports indicate that La/Al ratio 
in non-stoichiometric LAO films may be controlled dur- 
ing epitaxial growth. 13 ' 15 Atomically sharp interfaces are 
produced by molecular beam epitaxy (MBE) in which 
thermal energies of evaporated incident ions are low, 
about 0.1 eV, thus MBE avoids intermixing of cations at 
the interface. 5,11 However, the vaporization process used 
to facilitate transfer through the vapor LAO phase does 
not guarantee preservation of the target stoichiometry, 5 
which makes room for a possibility to control the film 
growth monolayer-by-monolayer. We note that for dif- 
ferent preparation methods of LAO/STO interfaces, e.g. 
pulsed laser deposition (PLD), other mechanisms may 
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give rise to conductivity. One of proposed mechanisms is 
formation of the high density of oxygen vacancies, which 
are generated in the STO substrate while depositing LAO 
thin film and can be responsible for increase of sheet car- 
riers density up to 5 x 10 17 cm~ 2 for PLD-grown n- 
LAO/STO interfaces if the sample is not annealed. 1 ' 16 
The insulating behavior of p-LAO/STO has been also 
ascribed to that the holes can be trapped by two elec- 
trons located at the oxygen vacancies created in the STO 
substrate. 17 

Yet another scenario for LAO/STO interface conduc- 
tivity that may take place in PLD-prepared structures 
is based on the suggestion that the La/Sr cation inter- 
mixing due to ion bombardment effect (inherent in PLD 
and post-growth treatment) may lead to the formation 
of one or two layers of metallic Lai-^Sr^TiOs. 5,18 ' 19 The 
thermodynamical stability for intermixed configurations 
has been recently reported. 5 ' 20 

In this paper, we aim to construct a clear picture 
of charge density redistribution both in stoichiomet- 
ric and non-stoichiometric interfaces of either type and 
LAO film thickness from 1 to 11 monolayers (0.5- 
5.5 u.c). The ab initio calculation methods employed 
are based on the density functional theory (DFT) us- 
ing a hybrid exchange-correlation functional. We con- 
trast stoichiometric/non-stoichiometric and p-type/n- 
type structures utilizing identical methods and compu- 
tation parameters. The B3PW functional 21 used in the 
CRYSTAL code 22 with atomic basis set (BS) contains a 
"hybrid" of the DFT exchange and correlation function- 
als with exact non-local Hartree-Fock (HF) exchange. 
For comparison, the selected set of interface configura- 
tions has been also modeled using the Perdew-Wang gen- 
eralized gradient approximation (PW91-GGA) density 
functional 23 ' 24 as implemented in the periodic plane- wave 
(PW) code vasp. 25 

We find that covalent effects in non-stoichiometric 
films are less pronounced than in stoichiometric ones 
and the structures are metallic in accordance with for- 
mal charges considerations. As Ti-0 bond strength ex- 
ceeds Sr-0 bond strength by ca. 120kJmol -1 (obtained 
considering formation enthalpies for respective oxides), 
in p-type IFs (where SrO monolayer is at the IF) we can 
expect covalent forces to be stronger than in n-type IFs. 
This, in its turn, leads to an approximately uniform free 
charge distribution through the film, whereas when cova- 
lent forces are weaker — as in n-type IFs — the free charge 
is forced to the edges of the LAO film (the surface and 
the IF) resulting in a bi-layered electron gas structure. 

Experimental works show that stoichiometric p-type 
interfaces exhibit no measurable conductivity, 11 but an- 
nealed stoichiometric n-type interfaces with LAO film 
thickness ^ 5 u.c. have free electron density in range 
1-3 x 10 13 cur 2 . 11 Similar densities (2-7 x 10 13 cm" 2 ) 
are obtained from first principle calculations 12-14 (cf. 
3.3 x 10 14 cm~ 2 predicted from electrostatic consider- 
ations). 

The paper is structured as follows. Section II describes 



the computational details of our calculations. The main 
part of the paper is formed by Sec. III. In Sec. Ill A we 
give an estimate of the thermodynamic stability and dis- 
cuss the electronic structure of ideal LaO- and AIO2- 
terminated LAO(OOl) surfaces. Section MB presents 
electronic charge distribution for n-LAO/STO and p- 
LAO/STO heterointerfaces and discusses their relation 
to the experimental and computational data available in 
the literature. Our conclusions are summarized in Sec- 
tion IV. 



II. COMPUTATIONAL DETAILS 

In this study LAO/STO heterointerfaces are modeled 
by means of two different methods: (i) linear combina- 
tion of atomic orbitals (LCAO) within the framework of 
hybrid density functional approach, and (ii) PW calcula- 
tions using the GGA density functional. 

To perform hybrid LCAO calculations, we used the 
periodic CRYSTAL code, 22 which employs Gaussian-type 
functions centered on atomic nuclei as the BSs for ex- 
pansion of the crystalline orbitals. The BSs used in this 
study were taken from the following sources: For Sr, Ti 
and O in the form of 311dlG, 411d311dG, and 8-411dlG, 
respectively, from Ref.26; for Al in the form of 8-621dlG 
from Ref.27; for La in the form of 311-31d3fl from crys- 
tal's homepage 22 (/-type polarization Gaussian function 
with the exponent a = 0.475 has been added according 
to prescription given in Ref.28). For Al and O all elec- 
trons are explicitly included. The inner core electrons of 
Sr and Ti are described by small-core Hay-Wadt effective 
pseudopotentials, 29 while the non-relativistic pseudopo- 
tential of Dolg et al. 30 was adopted for La. 

We employ the hybrid B3PW exchange-correlation 
functional 21 which accurately reproduces the basic bulk 
and surface properties of a number of AB0 3 perovskite 
materials. 26 > 31 ~ 33 The cutoff threshold parameters of 
CRYSTAL for Coulomb and exchange integrals evalua- 
tion (ITOL1-ITOL5) have been set to 7, 8, 7, 7, and 
14, respectively. Calculations were considered as con- 
verged only when the total energy obtained in the self- 
consistency procedure differed by less than 10~ 7 a.u. in 
two successive cycles. Effective charges on atoms as well 
as net bond populations have been calculated according 
to the Mulliken population analysis. 34-37 

As the second method the periodic total-energy code 
VASP 25 based on the use of a PW BS was applied. 
The cut-off energy has been chosen to be 520 eV. The 
non-local GGA exchange-correlation functional Perdew- 
Wang-91 (PW91) was employed. 23,24 Scalar relativis- 
tic projector augmented wave (PAW) pseudopoten- 
tials in our calculations contain 11 valence electrons 
(5s 2 5p 6 5d 1 6s 2 ) for La, 3 electrons (3s 2 3p 1 ) for Al, 10 
electrons (4s 2 4p 6 5s 2 ) for Sr, 12 electrons (3s 2 3p 6 3d 2 4s 2 ) 
for Ti, and 6 electrons (2s 2 2p 4 ) for O, respectively. Bader 
topological analysis 38 has been adopted to obtain net 
charges on atoms in VASP calculations. 
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TABLE I: Calculated equilibrium lattice constants (ao in 
A), atomic net charges (Qatom in e), cation-O bond popula- 
tions (-Pa/b-o m milli e), and band gaps (S in eV) of bulk 
LAO and STO in their high-symmetry Pm3m cubic phase. 
Shown are data obtained by means of both hybrid B3PW and 
standard GGA PW91 functionals. Negative bond population 
means atomic repulsion. Last two columns contain available 
experimental results for comparison. 
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5.51 
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In both VASP and CRYSTAL calculations the reciprocal 
space integration was performed by sampling the Bril- 
louin zone with the 8x8x1 Pack-Monkhorst mesh 39 
for all surface structures under consideration. For bulk 
computations we applied sampling with the 8x8x8 
Pack-Monkhorst mesh. Such samplings provide balanced 
summation in direct and reciprocal lattices. 

Taking into account that STO substrate at room tem- 
perature possesses perfect cubic structure, in our study 
we treat both LAO and STO in their high symmetry 
PmZm cubic phase. In fact, the bulk crystal structure 
of LAO, having space group i?3c (rhombohedral) with 
a = 5.364 A and c = 13.108 A at room temperature, 
can be represented by a pseudocubic unit cell with ao = 
3.790 A. At 821 K the structure of LAO transforms to 
become cubic with ao = 3.811 A. 40 Though the heteroin- 
terface assumes the transition between two intrinsically 
different crystal symmetries: Pm3m the substrate and 
R3c in the film, whereby thin films are expected to adapt 
to the substrate. 44 

Table I lists main bulk properties for both crystals. 
We note that the band gaps obtained by means of hy- 
brid B3PW computation scheme are in better agreement 
with experimentally observed results. Therefore in this 
paper we mainly discuss the results obtained by means 
of B3PW while results obtained using PW91 functional 
are published for comparative purposes in order to make 
our study consistent with earlier ab initio calculations 
performed basically on LDA- or GGA-DFT ground. 

Surface structures were modeled using a single slab 
model for LCAO calculations and a multi-slab model 
with vacuum gap of 20 A for PW calculation. To 
compensate the dipole moment arises at charged sur- 
faces, our slabs are symmetrically terminated. STO sub- 
strate contains 11 alternating (SrO)° and (Ti02)° atomic 
monolayers, while from 1 to 11 alternating (LaO) + and 
(A102)~ atomic monolayers were used for LAO film of 
the LAO/STO interface. Coordinates of all atoms in the 
LAO/STO heterointerfaces were allowed to relax. Due to 



TABLE II: Calculated deviations in surface monolayer net 
charge (AQ in e), and deviations of cation-O bond popula- 
tions (APa/b-o m milli e) in corresponding atomic monolayer 
relative to the bulk values (see Table I). Shown are data ob- 
tained by means of hybrid B3PW exchange-correlation func- 
tional. Surface monolayers are numbered beginning from the 
center of the slab (0 means the central monolayer of the sym- 
metrical slab unit cell). 
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symmetry constrains atomic displacements were allowed 
only along z-axis. Taking into account that the mismatch 
of ^2.5% between LAO and STO lattice constants arises 
during LAO epitaxial growth, in our modeling we have 
allowed relaxation of their joint lattice constant to mini- 
mize the strain effect. 



III. RESULTS AND DISCUSSION 

A. LAO(OOl) surfaces 

Before general discussion of LAO/STO interfaces stud- 
ied here, in this subsection we provide a comprehensive 
description of electronic and thermodynamic properties 
of both LaO- and A102-terminated pristine LAO (001) 
thin films. 



1. Electronic properties 

Pristine LAO(001) thin films were modeled using sym- 
metrical 9-monolayer slab model. Considering formal 
ionic charges, LAO(001) has alternating (LaO) + and 
(A102)~ surface monolayers and can be either LaO- 
or A102-terminated surface. Both LaO- and AIO2- 
terminations are studied. La/Al excess ratio is 1.25 and 
0.8 for LaO- and AlO 2 -terminated LAO(001) films, re- 
spectively. Monolayers in LAO(001) possess a net charge, 
the repeat slab unit cell has a non-zero dipole moment 
and therefore LAO(001) is type III polar surface accord- 
ing to Tasker's classification. 45 This means, that perfect 
and unreconstructed (1 x 1) LAO(001) surfaces consid- 
ered here can be stabilized by transferring of a half an 
electron (or hole) from the surface to the slab body that 
normally results in atomic and electronic reconfiguration 
at the surface. 

In Table II we list the changes in surface (LaO) + and 
(A102)~ monolayer net charges with respect to their 
bulk values (see Table I). Due to partly covalent nature 
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FIG. 1: (Color online) Projected density of states as calcu- 
lated by means of B3PW hybrid exchange-correlation func- 
tional: (a) LaO-terminated LAO(OOl), (b) A102-terminated 
LAO(OOl), (c) LAO bulk. TVB stands for the top of valence 
band. 

of La-0 and Al-0 bonds (positive Pa/b-o m Table I) 
net charges of La, Al, and O deviate from their formal 
ionic values of +3, +3, and —2, respectively. The La- 
O hybridization between La 5d and O 2p states lead to 
atomic charges of 2.43e, 2.07e, and — 1.50e for La, Al, 
and O, respectively. As a result, LaO and AIO2 mono- 
layers possess a bulk monolayer charge of ±0.93e instead 
of formal ionic ±le charge. According to the Table II 
surface monolayer of LaO-terminated LAO(OOl) attracts 
0.32 electrons, while other monolayers of the slab get the 
rest of 0.14 electrons to compensate the surface polarity. 
On the contrary, surface monolayer of AIO2 terminated 
LAO(001) solely receives 0.46 holes. Covalency of surface 
La-0 bond is only slightly increased (bond population in- 
creased only by 10 milli e), while calculated covalency of 
surface Al-0 bond is practically two times larger than in 
the bulk, that, to some extent, may compensate relatively 
modest surface relaxation of A102-terminated LAO(001) 
with respect to LaO-terminated one. 

Fig. 1 shows the density of states (DOS) projected onto 
all orbitals of La, Al, and O atoms of LAO bulk and both 
LaO- and A102-terminated LAO(001) surfaces as well. 
In case of LAO bulk (Fig. lc) the top of valence band is 
formed by O 2p orbitals, while the bottom of conduction 
band is formed mainly by La 5d states. La-0 hybridiza- 
tion is well pronounced. Calculated band gap of 5.51 
eV is in excellent agreement with its experimental value 
of 5.6 eV. 42 In case of LaO-terminated surface (Fig. la) 



gained excess of electrons shifts the Fermi level up to un- 
occupied level that gives raise to electron conductivity. 
In its turn the A10 2 -terminated surface (Fig. lb) expe- 
riences the lack of electrons that shifts Fermi level down 
to valence band and thus reveals the existence of hole 
conductivity. 

2. Thermodynamic stability 

The thermodynamic formalism adopted in the current 
study to estimate the stability of both LaO- and AIO2- 
terminated LAO(001) surfaces has been thoroughly de- 
scribed in Refs.46 and47 (see also references therein). 
The stable crystalline surface has to be in equilibrium 
with both LAO bulk and surrounding oxygen atmosphere 
assuming that an exchange of atoms between surface and 
environment is allowed. Therefore, the most stable sur- 
face has the lowest Gibbs free surface energy defined as 

MT,P) = j A [ ES t ah ~ N M Ek$ - (Nl, - iV A1 )A MLa 

-(N - 3N M )A» (T,p)}, (1) 

where t indicates the surface terminations, A the unit 
cell surface area, 2Vj the number of atoms of type i in 
the slab unit cell, Ef ab is the total energy of a slab 
with t surface terminations and is the LAO to- 

tal energy averaged per five-atom perovskite unit cell. 
A/Xj = [ii — -E^ulk' (i = La,Al) are deviations of chemical 
potentials for metal atoms from their energy in the bulk 
metals. For the oxygen atom such a deviation is consid- 
ered with respect to the energy of an oxygen atom in the 
ground triplet state of an O2 molecule Afio = /io~ ^-E-° 2 ■ 
Because pV term (V is unit cell volume) and the differ- 
ences in vibrational Gibbs free energy between the bulk 
solid and a corresponding slab is negligibly small, 48 we 
omit these two contributions. This permits replacing the 
Gibbs free energies in Eq. (1) and in the following for- 
mulas with the total energies obtained from ab initio cal- 
culations. 

In order to avoid the precipitation of relevant metals 
and oxides at LAO surface, as well as to prevent metal 
atoms to leave the sample the following conditions must 
be satisfied: 

> A/i La , > A,u A i, (2) 
eLmo 3 ~ E f Ah03 < 2A MLa + 3A Mo < < 2 o 3 , ( 3 ) 

where E' n is the formation energies of material n listed 
in Table III. 

We evaluate the oxygen chemical potential Afio(p,T) 
as a function of partial gas pressure and temperature us- 
ing the standard experimental thermodynamical tables 49 
as it was done in Refs.47,48. A^ (T,p 02 ) is the vari- 
ation of oxygen chemical potential due to temperature 
and pressure of the surrounding oxygen atmosphere. In 
addition to the experimental variation it contains a cor- 
rection term 5[ji = 0.03 eV, which compensates the dif- 
ference between the experimentally determined variation 
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TABLE III: Formation energies per formula unit used in 
analysis of surface stability. Experimental values are taken 
from Ref.49. 



Material 


E f , eV 


Exp. E f , eV 


La 2 3 


-17.52 


-18.64 


A1 2 3 


-16.68 


-17.37 


LaA10 3 


-17.68 





of the oxygen chemical potential and the reference state 
in current theoretical calculations (see Refs.50 and51 for 
a thorough discussion). 

Based on Eqs. 1, 2, and 3, the thermodynamic stabil- 
ity diagram is plotted in Figs. 2, showing the regions of 
stability of pristine LAO(OOl) surfaces with respect to 
precipitation of La 2 C>3 and AI2O3 oxides. Fig. 3 shows 
the thermodynamic stability diagram along the lines cor- 
responding to precipitation of La203 and AI2O3 oxides 
as a function of A/zo related to the temperature scale at 
an oxygen pressure typical during LAO/STO synthesis 
(P = 10 -6 mbar). To make such a diagram possible, ac- 
cording to prescription given in Ref.52 we replaced A/ZLa 
by 

AMLa=^(£L 2 o 3 -3AMo), (4) 

that corresponds to precipitation of La2 03 (lines 3 in 
Fig. 3) and by 

A/i La = £[ aA1G3 - \{E l M2Q3 - ^A Mo ), (5) 

that corresponds to precipitation of AI2O3 (lines 4 in 
Fig. 3). Formation energies for oxides are taken from the 
Table III. 

From the calculated thermodynamic stability diagrams 
we can predict that at ultra-high vacuum (UHV) con- 
ditions typical during PLD synthesis of LAO/STO in- 
terfaces and low temperatures (T < 550 K) the most 
stable is A102-tcrminatcd surface, while at elevated tem- 
peratures (T > 1100 K) stabilizes LaO-terminated sur- 
face. Between these temperatures both surface termina- 
tions may coexist. Further lowering of oxygen pressure 
shifts down these demarcated temperatures. This our 
prediction is in good qualitative agreement with time- 
of-flight scattering and recoiling spectrometry (TOF- 
SARS), atomic force microscopy (AFM), and photoelec- 
tron spectroscopy (PES) study performed by Rabalais 
and co-workers. 53 ' 54 They found that at temperatures 
less than 423 K, the surface is exclusively terminated by 
an Al-0 layer, while at temperatures above 523 K the 
surface is exclusively terminated by a La-0 layer. Be- 
tween 423 K and 523 K surface stoichiometry changed 
from AlOz to LaO^ and thus mixed terminations were 
proposed. Moreover this change was found to be fully 
reversible. Rabalais and co-workers suggested that the 
surface termination change was caused by the forma- 
tion of surface oxygen vacancies at high temperature, 



which drives the migration of the La atom to the sur- 
face and the Al atom into the bulk. More recent ex- 
perimental study based on X-ray crystal truncation rod 
(CRT) analysis 55 demonstrates that LAO(OOl) possesses 
Al-terminatcd structure at both room and high (670 K) 
temperatures with no evidence for the reversal of sur- 
face termination or for the formation of surface oxygen 
vacancy. Authors of Ref.55 explain the observation of La- 
rich termination in ion-scattering experiments 53 ' 54 by the 
effect of the increasing access to the lanthanum atom be- 
cause of considerable surface oxygen relaxation that leads 
to a significant enhancement of the lanthanum atom sig- 
nature. On the other hand, Marx and co-workers have 
observed the La-terminated LAO (001) with stoichiome- 
try of (VLa^s) -05 , where V is the lanthanum cation 
vacancy, i.e., each surface La is coordinate to four sur- 
face oxygens and four oxygens in the subsurface layer. 56 
Therefore one may conclude that the experimental analy- 
ses have been performed at various conditions and report 
either LaO- and A10 2 -terminated LAO (001) or mixture 
of them, so it is not clear if surfaces reached thermody- 
namic equilibrium or not. 

Ab initio thermodynamical stability diagrams pre- 
viously calculated for LAO(001) shows that LaO- 
terminated surface is more stable with respect to AIO2- 
terminated one 57 and LaO-terminated surface contain- 
ing oxygen vacancy is more stable than oxygen defi- 
cient A10 2 -terminated LAO(001) as well. 58 Mixed sur- 
faces with LaO- and A102-terminations were not pre- 
dicted. In fact, our thermodynamic analysis does not 
support this prediction. From our point of view the main 
reason for such a discrepancy may be the different com- 
putational approach, DFT within local density approxi- 
mation, used by Authors of Refs.58,57. 



B. LAO/STO heterointerfaces 

1. Charge redistribution and electronic properties 

Calculations of electronic properties of the 
LAO/STO(001) heterointerfaces were carried out 
using the symmetrically terminated slab model. The 
STO(001) substrate consisted of 11 atomic monolayers 
and could be terminated with either (Ti02) monolayer 
in n-type heterostructures or with (SrO) monolayer in 
p-type heterostructures. Then monolayer- by-monolayer 
epitaxial growth was modeled adding a pair of respective 
monolayers of LAO(001) symmetrically to both sides of a 
substrate slab until deposited LAO(001) thin film reach 
thickness of up to 11 monolayers. In such way we con- 
struct 22 heterostructures of both types and of different 
LAO film thickness to model. Note that 11-monolayer 
thick substrate and 20 A thick vacuum gap used for 
vasp GGA calculations is enough to avoid undesirable 
interaction of neighboring surfaces/interfaces and allows 
us to reach the equilibrium charge density redistribution 
in heterointerfaces under study. Due to the restrictions 
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FIG. 2: (Color online) Thermodynamic stability diagram as a function of O and La chemical potentials built for both LaO- 
and A102-terminated LAO(001) surfaces. Diagram contains precipitation conditions for both La and Mn metals, as well as for 
their trivalent oxides (La203 and AI2O3). Stable region is shown as shaded area between La203 and AI2O3 precipitation lines. 
The numbers from 1 to 4 in the circles indicate segregation lines for precipitation of: 1. Al, 2. La, 3. La2C>3, 4. AI2O3. The 
right side shows a family of oxygen chemical potentials under different conditions. The label m indicates the O2 gas partial 
pressure: 10 m mbar. Red (gray) line corresponds to oxygen partial pressure p = 0.2po as in the ambient atmosphere. Point 
A stands for room temperature and ambient oxygen pressure, point B stands for typical temperature and pressure during 
LAO/STO(001) synthesis. 
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FIG. 3: (Color online) The thermodynamic stability di- 
agram calculated along the La203 and AI2O3 precipitation 
lines (number 3 and 4 in the circles, respectively) with A/iLa 
defined according to Eqs. 4 and 5. The dependence on 
the oxygen chemical potential is converted to the appropri- 
ate temperature scale at an oxygen pressure typical during 
LAO/STO(001) synthesis (P = 10" 6 mbar). The interval be- 
tween points C and D correspond to temperature range where 
both LaO- and A102-terminated LAO(001) surfaces are stable 
and may coexist. 



by imposed symmetry, in our calculations atomic posi- 
tions of all the heterointerfaces under study were relaxed 
along the z axis. 

If we consider atomic displacements, we can see that 
cations and anions in LAO monolayers have considerably 
different displacements, thus electric dipole moment ap- 
pears and accumulates within the thin film. Stoichiomet- 
ric heterointerfaces have greater displacement differences 



between anions and cations than non-stoichiometric ones 
in LAO monolayers, while the situation is diametrically 
opposite for the STO monolayers. As we shall see fur- 
ther, the dipole moment creates an electric field, and its 
potential strongly correlates with the distortion of the 
band edges (so-called polar distortion), which then gives 
rise to the conductivity in stoichiometric LAO/STO(001) 
heterointerfaces of n-type. 

To predict the charge redistribution in heterointerfaces 
we calculated the changes of net atomic Mulliken charges 
in comparison with the bulk phase of the LAO and STO 
parent materials. These charge deviations are shown in 
Fig. 4a-d for LAO/STO(001) heterointerfaces of n- and 
p- type. From these one can clearly see, that deviation of 
charges are relatively small in the inner monolayers of the 
LAO film in n-type LAO/STO(001), not exceeding 0.03e, 
whereas the same layers in the p-type LAO/STO(001) 
show quite large charge deviations ±(0.35-0.40)e from 
the parent bulk, and these are negative for AIO2 mono- 
layers and positive for LaO monolayers. 

In both n- and p-type interfaces charges on the sub- 
strate monolayers did not vary substantially. For sto- 
ichiometric n-type and non-stoichiometric p-type inter- 
faces these are about ±0.04e for Ti02 and SrO, respec- 
tively. On the other hand, stoichiometric p-type inter- 
faces show a small positive deviation of Ti02 monolayer 
charges (ca. O.Ole) and about ten times bigger negative 
charge deviation for SrO monolayers. Charge shifts in 
the substrates of stoichiometric n-type structures are all 
negative, and SrO shifts (ca. 0.04e) are smaller than Ti02 
shifts of ca. 0.06e. 

Most significant deviations in atomic charges of n- 
type structures are located in the top-most monolayer — 
+0.2e for stoichiometric structures and — 0.25e for non- 
stoichiometric ones — due to the surface effects and thus 



7 



o. 
o 




-0 
-0 

'oT 

w -0. 

a 

< o. 





-0 
-0 
-0 
-0. 



(a) 



(b) 



(c) 



(d) 



0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 

0.4 
0.3 
0.2 
0.1 


-0.1 
-0.2 
-0.3 
-0.4 



o <n <t r- pn in c^r- f o w n ^ >c oo o*" 
Atomic Plane Atomic Plane 

FIG. 4: (Color online) Calculated deviations of Mulliken ef- 
fective charge densities (APq) in AO- and B02-monolayers 
of (a, b) n-LAO/STO(001) and (c, d) p-LAO/STO(001) het- 
erostructures with respect to charge densities in AO- and 
B02-monolayers of STO and LAO bulk, correspondingly. 
Calculations are performed using B3PW hybrid exchange- 
correlation functional. The x axis shows the atomic monolay- 
ers from which atoms are originated. STO and LAO monolay- 
ers are numbered starting from the center of slab (0 means the 
central monolayer of the symmetrical slab unit cell). Mono- 
layers (planes) are numbered separately for STO(001) sub- 
strate and for LAO(001) nanofilm. Panels (a, c) show charge 
density deviation for A?lao = 10, while panels (b, d) — for 
JVlao = 11. 



compensate the "polar catastrophe" as proposed from a 
pure ionic model. 59 In p-type structures charge shifts in 
the surface layers are less pronounced than in the inner 
layers of the film and are +0.05e and — 0.27e for LaO- 
and A102-terminated structures, respectively. 

Here charge redistribution only in the thickest struc- 
tures investigated is shown. Respective graphs for thin- 
ner structures can be found in Re fs. 60,61. 

Another way to look at the problem of charge redistri- 
bution is to calculate, what happens with the electronic 
charge density in the heterostructures, compared to the 
isolated LAO and STO slab parts. Charge density redis- 
tribution is defined as the electronic density in the het- 
erointerface minus the sum of electron densities in sep- 
arately isolated STO(001) substrate and LAO(001) thin 
film slabs and is depicted in Fig. 5 for both n- and p-type 
LAO/STO(001) interfaces. 



These plots show us that the most significant distor- 
tions occur at the interface due to the compensation of 
the surface effects of the slabs. They also show that the 
electronic structure of the substrate of non-stoichiometric 
heterostructures is distorted stronger than that of stoi- 
chiometric ones. The situation in the thin films is oppo- 
site. This fact correlates with the argument in the section 
on atomic structure. 

More illustrative property to consider is the polariza- 
tion of all four of stoichiometric and non-stoichiometric 
n- and p-type heterointerfaces, which was already briefly 
introduced. It allows us to explain certain phenomena, 
such as the polar distortion, as well as to provide a mech- 
anism for a partial compensation of the "polar catastro- 
phe". 

Let us assume that each one of considered interfaces 
possesses no net charge, thus it can be divided into mul- 
tiple neutral slabs normal to z, in which net charge is 
also zero and average polarization of such slabs can be 
calculated. Charge density function that should be used 
in the calculations, is estimated as if the charge of each 
atom A is uniformly distributed over the plane z = za, 
reducing the task to one dimension. Thus the projection 
of polarization vector on z-axis can be calculated as 



Pi 



Ea z aQa 
Az • 



(6) 



Qa is the charge on atom A, Az is the thickness of the 
neutral allocated slab, to which the atom A belongs and 
summation is performed over all the atoms in the z-th 
neutral slab. 

In order to divide the interface in neutral slabs, it 
sometimes is necessary to split one monolayer's charge: 
One part of it compensates the remaining charge of the 
previous slab and the remainder goes to the next one. 

The resulting polarization function P(z) is averaged 
using the moving average function, and the results for n- 
and p-type interfaces are shown together with the ener- 
gies of band boundaries Stvb and -Ebcb and the poten- 
tial due to intrinsic electrostatic field V in Fig. 6. Here 
one can see, that LAO films of stoichiometric interfaces 
are strongly polarized, giving rise to the polar distor- 
tion of band edges. On the other hand, there is rather 
weak LAO polarization in the non-stoichiometric inter- 
faces meaning a weak polar distortion as is observed. The 
substrate is polarized more in non-stoichiometric case, 
which corresponds to Figs. 5c, d and Figs. 5g,h. AIO2- 
terminated structures possess substantial polarization in 
the top-most monolayer. Top-most layer's polarization 
of LaO-terminated structures, on the other hand, is neg- 
ligible. The interface monolayers of n-type structures are 
substantially polarized. 

Electronic properties in a more experimentally mea- 
surable way can be represented as band gaps for insu- 
lating structures or as the concentration of charge car- 
riers for conductors. These data obtained with CRYS- 
TAL and vasp are represented in Tables IV and V for n- 
and p-type structures, respectively. Firstly, one can see 




FIG. 5: (Color online) Difference electron charge density maps calculated for (a-d) n-LAO/STO(001) and (e-h) p- 
LAO/STO(001) heterostructures: (a, e) (110) cross-section for Alao = 10, (b, f) (100) cross-section for Alao = 10, (c, 
g) (110) cross-section for Alao = H, (d, h) (100) cross-section for -/Vlao = 11. Red (dark gray), blue (light gray) and gray 
isolines describe positive, negative and zero values of the difference charge density, respectively. Isodensity curves are drawn 
from -0.025 to +0.025 e A~ 3 with an increment of 0.0005 e A" 3 . Right-side bar shows the atomic monolayers from which atoms 
are originated. Calculations are performed using B3PW hybrid exchange-correlation functional. STO and LAO monolayers 
are numbered beginning from the center of slab (0 means the central monolayer of the symmetrical slab unit cell). Monolayers 
(planes) are numbered separately for STO(001) substrate and for LAO(001) nanofilm. 



that all the non-stoichiomctric interfaces arc conducting 
and free charge concentration is roughly equal within a 
type and does not depend on the LAO film thickness. 
p-Type structures possess greater carrier density than n- 
type structures, though experiments never showed con- 
ductive behavior in the former. 

For stoichiometric structures insulating behavior is the 
default one. The thickness of the band gap decreases with 
the thickness of the LAO film both for n- and p-type 
structures. This eventually leads to the closing of the 
gap for the n-typc interfaces with A^lao ^ 10 monolayers, 
which is in a good accordance with experimental works. 11 
The gap-diminishing tendency is less pronounced for the 
p-type structures and thus they are not found conducting 
at any thickness within this study. 



The results obtained with VASP are given for qualita- 
tive comparison. They showed out to be in accordance 
with crystal results, but due to the specifics of the non- 
hybrid functional band gaps and free charge concentra- 
tions are far too small. Taking into account that the 
largest difference between calculated using CRYSTAL code 
and experimentally observed band gap of bulk materials 
is 0.39 eV (see Table I) we note that our crystal calcu- 
lations give plausible results comparing to experimental 
data. 

The total band gap described above gives us some valu- 
able data on conducting-insulating behavior of the inter- 
faces of different types. Nevertheless, it does not give 
us much information about the origin of conductivity. 
Thus it is more worthy to look at the positions of the 
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FIG. 6: (Color online) Polarization (as calculated using 
Eq. (6)), band edges and electrostatic potential of (a, b) 
n-LAO/STO(001) and (c, d) p-LAO/STO(001) heterostruc- 
tures with (a, c) Alao = 10 and (b, d) Alao = H LAO 
monolayers. Zero at the energy scale corresponds to the Fermi 
level. Distances are measured from the central monolayer of 
the symmetrical slab unit cell. TVB stands for the top of va- 
lence band, BCB stands for the bottom of conduction band. 



TABLE V: The same as Table IV, but for p-LAO/STO(001) 
heterostructures. 







B3PW (crystal) 


PW91 (vasp) 


atLAO 
Jv tot 


Term, m-layer 


8 




5 n s 


1 


AlOa 




6.65 


- 


2 


LaO 


4.00 




1.60 


3 


AlOa 




7.27 


- 


4 


LaO 


4.05 




1.69 


5 


A10 2 




9.08 




6 


LaO 


4.05 




1.51 


7 


AlOa 




7.90 




8 


LaO 


3.80 




0.48 


9 


A10 2 




6.97 




10 


LaO 


2.92 




0.25 


11 


AlOa 




10.2 


0.12 



than investigated. Non-stoichiometric interfaces show 
little or no polar distortion, but it is not necessary for 
the appearance of the conductivity, because such struc- 
tures contain non-stoichiometric LAO films, which are 
already conducting on their own. Our prediction on con- 
ductivity of non-stoichiometric LaO-terminated n-type 
LAO/STO(001) interface is in agreement with a recent 
theoretical study performed by Pavlenko and Kopp (See 
Ref.62) in which they show that LaO-terminated n-type 
LAO/STO(001) interface is metallic. 



TABLE IV: Band gaps (5 in eV) or sheet carrier density (n s 
in 10 14 cm' 2 ) of n-LAO/STO heterointerfaces as calculated 
by means of hybrid B3PW and PW91 exchange-correlation 
functionals. stands for the total number of LAO(001) 

monolayers deposited atop STO(001) substrate. 



B3PW (crystal) PW91 (vasp) 



/V LAO 

iv tot 


Term, m-layer 


5 


n s 


5 


n„ 


1 


LaO 




6.04 






2 


AlOa 


3.65 




1.41 




3 


LaO 




6.07 






4 


AlOa 


2.91 




1.03 




5 


LaO 




5.91 






6 


AIO2 


1.96 




0.40 




7 


LaO 




6.20 






8 


AIO2 


1.07 




0.03 




9 


LaO 




6.27 






10 


AlOa 




1.56 




0.16 


11 


LaO 




6.13 




0.54 



band edges in energy scale separately for each mono- 
layer. Such a decomposition is depicted in aforemen- 
tioned Figures 6a, b and 6c, d for n- and p-type structures, 
respectively. From these plots one can see, that band 
edges for stoichiometric interfaces are distorted, besides 
such a distortion leads to n-type conductivity in n-type 
structures thick enough and might hypothetically lead 
to the p-type conductivity in thicker p-type structures 



IV. SUMMARY AND CONCLUDING 
REMARKS 

We have performed large-scale first-principles calcu- 
lations on a number of both stoichiometric and non- 
stoichiometric LAO/STO(001) heterostructures. Two 
different ab initio approaches have been applied: LCAO 
with hybrid B3PW and PW with PW91 exchange- 
correlation functionals within DFT. Consistently within 
both approaches we predict that there exists a distortion 
in energies of band edges for stoichiometric structures 
which eventually leads to the appearance of the conduc- 
tivity at a critical thickness in n-type interfaces or to the 
reduction of the band gap for p-type interfaces. Non- 
stoichiometric interfaces were found to be conducting in- 
dependently of the LAO film thickness and possessing 
little or no distortion of band edges. The conductivity 
appears due to the non-stoichiometry of the thin film 
which is a conductor on its own, as we demonstrate by a 
separate analysis of an isolated film. 

The degree of distortion of the band edges agree well 
with the estimates of the internal electric field generated 
by changes in the atomic charges and the geometric relax- 
ation of the atomic structure. We confirm these factors 
as the ones responsible for the rise of conductivity in sto- 
ichiometric n-type heterostructures. Calculated concen- 
tration of the free charge in the interfaces roughly agrees 
with the experimental data, being somewhat underesti- 
mated. 
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For non-stoichiometric n-type interfaces electron gas 
structure is monolayered with uniform distribution over 
both the film and the substrate, while for p-type inter- 
faces it is bilayered with one part of free charge carriers 
located on 3d orbitals of Ti at the IF, while the other is 
located on La orbitals at the surface. The total calculated 
n«6x 10 14 cm -2 well accords with that, predicted from 
electrostatic assumptions tt-es = 1/a 2 = 6.6 x 10 14 cm~ 2 . 
Of that the IF gas layer gets nrp ~ 1.3 x 10 14 cm~ 2 and 
the surface gas layer gets ns ~ 4.7 x 10 14 cm~ 2 . 

Thermodynamic analysis that we have performed for 
the pristine LAO(OOl) surface reveals that its both LaO- 
and A102-terminations may co-exist at temperatures 
above 550 K. If LAO/STO(001) heterointerface is cov- 
ered by LaO monolayer, charge compensation mechanism 
of deposited polar non-stoichiometric LAO film leads to 
the tendency of Ti 3+ formation at the interface (see Fig. 
4). To some extend it may explain the unexpected ob- 
servation of Ti 3+ photoemission spectroscopy peak from 
n-type LAO/STO interfaces grown at 873 K. 63 

In general, we conclude that one should not disre- 



gard the stoichiometry aspect when considering ways 
to make the LAO/STO interfaces conducting as non- 
stoichiometric interfaces possess unique quasi-2D elec- 
tron gas structure that gives an overall 2 times greater 
free charge carrier density in comparison with stoi- 
chiometric interfaces. For stoichiometric n-type struc- 
tures covalent and electrostatic forces' interplay leads to 
metal insulator transition at critical film thickness, but 
for non-stoichiometric — to formation of bilayered (n-type 
IFs) or monolayered (p-type IFs) quasi-2D electron gas. 
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FIG. 7: (Color online) Calculated deviations of Mulliken effective charge densities (APq) in AO- and B02-monolayers 
of n-LAO/STO(001) heterostructures with respect to charge densities in AO- and B02-monolayers of STO and LAO bulk, 
correspondingly. Calculations are performed using B3PW hybrid exchange-correlation functional. The x axis shows the atomic 
monolayers from which atoms are originated. STO and LAO monolayers are numbered starting from the center of slab (0 
means the central monolayer of the symmetrical slab unit cell). Monolayers (planes) are numbered separately for STO(001) 
substrate and for LAO(001) nanofilm. Number of monolayers increases from 1 to 11 for panels (a) to (k). 
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FIG. 8: (Color online) Calculated deviations of Mulliken effective charge densities (APq) in AO- and B02-monolayers 
of p-LAO/STO(001) heterostructures with respect to charge densities in AO- and B02-monolayers of STO and LAO bulk, 
correspondingly. Calculations are performed using B3PW hybrid exchange-correlation functional. The x axis shows the atomic 
monolayers from which atoms are originated. STO and LAO monolayers are numbered starting from the center of slab (0 
means the central monolayer of the symmetrical slab unit cell). Monolayers (planes) are numbered separately for STO(OOl) 
substrate and for LAO(OOl) nanofilm. Number of monolayers increases from 1 to 11 for panels (a) to (k). 



